Phase-reversed structures in superlattice of nonlinear materials 
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We present detailed description of so-called phase-reversed structures that are characterized by 
two grating wave vectors allowing simultaneously phase- match two parametric three- wave processes. 
The novelty is that the structure is realized as a definite assembly of nonlinear segments leading to 
detailed description of cascaded three-photon processes with the parameters of realistic structured 
nonlinear materials of finite length. We apply these results for analysis of the quasi-phase-matching 
in production of both photon triplet and four-photon states in cascaded down-conversion. The 
received results are matched with the experimental data. 
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I. INTRODUCTION 

Cascading processes in optical materials that involve 
several different second-order nonlinear interactions pro- 
vide an efficient way to lower the critical power of optical 
parametric devices. The concept of multistep cascading 
brings new ideas into this field, leading to the possibility 
of an enhanced nonlinearity- induced phase shift, the si- 
multaneous generation of higher-order harmonics, multi- 
color parametric interaction, difference- frequency genera- 
tion, (see, [l[ and references therein) . Multistep paramet- 
ric interactions and multistep cascading are characterized 
by at least two different nearly phase-matching or quasi- 
phase-matching (QPM) parameters. QPM is an impor- 
tant technique in nonlinear optics and is widespread used 
for various applications 0, 0] ■ The width of the phase- 
matching curve depends on the type of the crystal, its 
length, and the phase-matching method. QPM not only 
makes efficient frequency conversion possible, but also en- 
ables diverse applications such as beam and pulse shap- 
ing, multi-harmonic generation, high harmonic genera- 
tion and all-optical processing. Recently, it has been 
demonstrated that the quantum-computational gates, 
multimode quantum interference of photons and prepara- 
tion of multi-photon entangled states also can be realized 
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FIG. 1: Scheme of periodically poled nonlinear crystal 
(PPNC) which involves nonlinear domains of length / with 
X+ > and %- < with phase-reversal grating (with period 
Aph) upon a uniform QPM grating (with period d). 
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in QPM structures [3H6|. 

One of the non-uniform QPM structures is so-called 
phase-reversed structures [7J, which has been studied 
both theoretically and experimentally and have been 
employed for the wavelength conversion, and for the 
cascaded single-crystal third harmonic generation pro- 
cess @, Q. Several proposals for increasing the phase- 
matching region have also been suggested, including the 
use of the phase-reversed QPM structures. These struc- 
tures are based on the use of a phase-reversed sequence 
of many equivalent uniform short QPM domains, par- 
ticularly, periodically poled domains, in which neighbor- 
ing nonlinear segments have the opposite signs of the 
second-order susceptibility. The phase-reversed struc- 
tures arranged in such a way that at the place of the 
joint two ends of the domains segments have the same 
sign of the x coefficient (see Fig. I). On the whole, 
this structure is characterized by two grating wave vec- 
tors allowing phase-match two processes. In the standard 
approach, the one-dimensional x^ 2 \ z ) function leading 
to the phase-reversed QPM structure can be described 
by the infinite Fourier expansion on two periods. Such 
phenomenological approach qualitatively explains the si- 
multaneous QPM of two parametric nonlinear processes, 
however, has not included the parameters of realistic 
structured nonlinear materials. These dates often should 
be taken into to consideration if we are interest in more 
detailed description of cascaded parametric processes. 

In this paper, we present more detailed description of 
phase-reversed structures in another way as the definite 
assembly of nonlinear segments. Then, we arrange this 
structure for two cascaded parametric devices proposed 
for generation of three-photon and four-photon states. 
In this way, we consider cascaded collinear parametric 
processes in one dimensional media. The paper is ar- 
ranged as follows. In Sec. II we consider an effective 
nonlinear coefficient G(Afc) of three-wave interaction in 
one-dimensional phase-reversed media. This quantity is 
defined as the Fourier transform of the second order sus- 
ceptibility x^ 2 ' l ( z ) hi term of wave vectors mismatch Ak 
of three- wave interaction caused by dispersion in an opti- 
cal material. In Sec. Ill we briefly consider application of 
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these results for obtaining QPM in generation of multi- 
photonic states. 



II. PHASE-REVERSED STRUCTURES OF 
FINITE LENGHT 



where Afc is the phase-matching parameter and L is a 
length of medium. In the phenomenological approach 
the x^ 2 \ z ) function leading to the phase-reversed QPM 
structure can be described by the following Fourier ex- 
pansion 



We start with consideration of effective coupling con- 
stant G(Afc) that usually describe three- wave interaction 
in second-order nonlinear media with the susceptibility 
X^ 2 \ z ) m the following form 



G(Afc) 
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that is the product of two periodic functions. 
go = 0, 5„^ = 2/7rn, g m ^ = 2 firm and G„ 
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= -r^-m. In this case the effective interaction con- 
stant reads as 
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Considering the phase-reversed QPM structure as the 
definite assembly of nonlinear segments we turn to the 
general formula obtained for the effective coupling con- 
stant in multilayered media [loj 



k™(z,Lj) — ^-n n (z,ujj), n n {uj) is the refractive index of 
the medium at the given frequency. 

In the result, the effective nonlinear coefficient G(Afc) 
for the scheme depicted in the Fig.l can be calculated in 
the following form 



G (Afc) = J2 l m X m e-^ + ^) Sl nc f^Hl) , 



(4) 



where Ak n = kfi — fc™ — fcj is the phase-matching func- 
tion of three photon process in the each n-th domain, 



G(Ak) = L Xo e- ia ^Y M ,N(Ak). 



(5) 



We assume that phase-matching parameters are the same 
for all layers Ak n = Ak, as well as x+ = \X- \ = Xo- 

We introduce the QPM function YM,N(Ak) which is 
the product of the standard phase-matching spectral 
function for a single segment, i.e. s 
assemble function 



'^), and the 
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Here: iV-is the number of domains in each block, M-is 
the number of blocks, a(Afc) = ^[(Afc - G) + [M - 
l)(Afc — F)], l-is the length of each domain, G = 2w/d, 
F = 2tt/A p , 1 , d = 21, A ph = 2NI, L = MNl. It is easy to 



check that the function YM,N{Ak) has an alternative rep- 
resentation, as the product of expansions on the number 
of domains and blocks 
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where tp — |(Afc — G). As we see, the above expressions 
(HJ), fl7j derived on the base of the method of super latices 
elaborated on the multilayered media differs from the for- 
mula ([3]) , which is obtained in the standard, phenomeno- 
logical approach. In contrast to the expression (O, the 
results ©,([71) depend on the number of the domains and 
blocks and also describe the case of small numbers of the 
segments. It is remarkable, that in this case the QPM 
function displays specific interference effects. Neverthe- 
less, as shows our analysis some consequences of both ([3j) 
and results, for N 3> 1, M » 1, coincide qualita- 
tively. 

Note, that single-phase-matched process of three- 
photon interaction in multilayered media is also described 
by an assembly function. For PPLN structure (see, one 
of the block in Fig.l containing N layers) the phase- 
matching function is calculated as 
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FIG. 2: Plot of the function Ym,n(x) for even numbers of 
both domains and blocks; N = 22 and M = 8. The inset 
shows the twin narrow peaks of opposite amplitudes. 
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where Afc is the wave vectors mismatch of three-wave 
interaction. This function exhibits a structure of single 
narrow peaks as discussed in (Tlj . 

The obtained QPM function Y"M,iv(Afc) allows simul- 
taneously phase-match two parametric processes which 
is achieved by using two grating vectors. In order to il- 
lustrate the possibility of the double phase-matching in 
such superlattice structure we present a detailed analysis 
of the function Ym,n{x) in dependence on the parame- 
ter x — i n the graphical form (see, Figs. [2J [3]). The 
assembly function exhibits the structure of peaks that 
are separated by n. Unlike to the assembly function 
for single-phase-matched process, Yjv(Afc) function has 
more complex structure, particularly, each group of the 
extremes consists of thin structure of twin-narrow peaks, 
separated by oscillating low maxima. The sine (^^) 
multiplier causes decreasing the amplitude of the con- 
secutive peaks. The behavior of Ym,n(%) function is dif- 
ferent for even and odd values of M. As depicted in 
Fig. [5] for even values of M the twin peaks have the op- 
posite signs, while for odd values of M the both narrow 
peaks have the same sign (see Fig. [3j. The existence of 
twin peaks and oscillation between them are the main 
peculiarities of this structure. These twin narrow peaks 
can be chosen to phase-match two parametric processes 
involved into the multistep-cascading. 

Here we represent Ym ! jv(AA;) function for the small 
number of domains to exhibit direct correlation between 
number of domains and the width of phase-matching 
function (see, Fig. [4]). It is easy to notice that with the de- 
crease of number of domains and blocks the narrow peaks 
in the pase-matching function becomes wider, as well as 
the number of oscillation between twin peaks decreases. 
Besides this, the positions of the maxima of the phase- 
matching function becomes slightly shifted, although the 
periodicity of the Ym,n{x) function is remained the same. 




FIG. 3: Plot of the function Ym,n(x) for even numbers of 
N = 22 and odd numbers M = 9. The inset shows the twin 
narrow peaks of the same amplitudes. 



III. APPLICATIONS: QPM FOR PRODUCTION 
OF MULIPHOTON STATES 



In this section, we apply the above results analyzing 
the QPM in generation of both photon triplet and four- 
photon states in cascaded down-conversion. Recently, 
production of joint multiphoton states, particularly, pho- 
ton triplet have attracted a great interest in probing the 
foundations of quantum theory and for various applica- 
tions in quantum information technologies (see, |12Hl5l | 
and the references therein). Depending on the configu- 
ration of the experiment and the designs of the crystal, 
the multiphoton states can be entangled in any num- 
ber of variables: time, frequency, direction of propaga- 
tion, and polarization. Studies carried out with miltipho- 
ton state beams range from the examination of quantum 
paradoxes, to applications in optical measurement s, sp ec- 
troscopy, imaging, and quantum information (see, [l3lll6j 
and references there). In this way, the direct generation 
of photon triplets using cascaded photon-pairs has been 
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FIG. 4: Plot of the function Ym,n(x) for the small numbers 
of domains and blocks, namely N = 6 and M — 2. The insert 
shows the broadening of the peaks. 



demonstrated in periodically poled lithium niobate crys- 
tals [l4j |. Most recently, joint quantum states of three- 
photons with arbitrary spectral characteristics have been 
studied on the base of optical superlattises [l5[ . In con- 
tinuous of this paper, we analyze the cascaded generation 
of photon triplets in phase-reversed superlattice of x^ 
materials. 

We consider collinear, one-dimensional configurations 
focusing on consideration of cascaded parametric down- 
conversion (PDC), which generates three-photon states 
in the presence of an optical cavity. In the scheme con- 
sidered the pump field at the frequency ioq converts to 
the subharmonics at the central frequencies to± = ^ 
and L02 = ^f 1 throughout two cascaded processes: to a — >• 
lji + u>2 and u>2 — > to\ + u)\. Below we consider the am- 
plitude of three-photon generation through the cascaded 
processes assuming the frequencies w\ and u>2 as variable 
quantities distributed around the central values ^ and 
respectively. These processes are characterized by 
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the wave- vector mismatches: Afci = k(t0o) — k(uii) — k(ui2) 
and Ak2 = k(uj2) — — The model Hamilto- 

nian for the system, in the rotating-wave approximation, 
is given by 

H mt = ih (E e- lulot a+ - E*e lu)ot a ) + 
ihC (a a+b + - a+ab) + ih£ (a +2 b - a 2 b+) . (9) 

Here, £ and £ are the nonlinear coupling constants that 
are the Fourier transformations of the second-order non- 
linear susceptibility x^H 2 ) ( see ; Eq-CE})- The cavity 
modes are described by discreet creation and annihila- 
tion operators a + ,b + and a, b as Eq describes the am- 
plitude of the driving fields. For the case of the phase- 
reversed superlattice structure the constants of the pro- 
cesses ujq — > lu± + u>2 and u>2 — > &i + are calculated 
as 



where x\ = lAki/2, X2 = lAk2/2. In general, the ampli- 
tude of three-photons at the frequencies io\ , cj 1; to 1 gen- 
erated in cascaded down-conversion has been calculated 
in [l^]. This quantity is proportional to the factor -£?oC£- 
For the case of phase-reverse structure we can calculate 
three-photon amplitude in the following form 



$(wi, WijWi) = f(ui,oJi)h(Aki, Ak 2 ), 



(11) 



where 



C~ie- iQ ( Afel )F M ,iv(ari), 
£ ~ Le' la ^Y M , N (x 2 ), 



(10) 



h(Ak 1 ,Ak 2 ) = {MNfY M ^{k 1 )Y M ^(k2) (12) 

is the two-dimension assemble function that describes 
joint phase-matching, while the function f (101,1*12) — 
Cl 2 xoE (ui )e~ ia(Akl) e' la( - Ak2 \ depends on the param- 
eters of three-wave interaction, the phase multipliers 
a(Aki) = ^[(Aki-G)+(M-l)(Aki-F)} in the formula 
§5§ and C is the normalization factor of three-photon 
amplitude. Three-photon down-conversion is controlled 
by energy conservation between the pump and daughter 
photons uiq = u>i + lo 1 + to 1 . 

Analyzing the graphics of Ym,n(x) (see Fig. [3]) in appli- 
cation to cascaded scheme of down-conversion we could 
arrange two twin peaks of the phase-matching function 
with the wave- vector mismatches of cascading processes. 
We demonstrate it using the experimental dates for the 
cascaded down-conversion [lj|- For the typical value 
of pump wave length Al = 0.53/jto and LiTaC>3 crys- 
tal mismatch parameters from experimental results are 
Afci = 0.32/mi- 1 and Ak 2 = 0.87/jm" 1 . We check 
that the simultaneous phase-matching of two processes 
is really realized on the formula (|6|) for the parameters: 
I = 10.25/jm, N = 22, M = 9. In this case, the Afci 
match with the second peak of the first twin-peaks on 
the positive branch of the axes and the Ak2 matches 
with the first peak of the second twin-peaks, which are 
situated next to the first twin-peaks and separated from 
them by the period. We got perfect correspondence of 
theoretical and experimental results. 

The dependence of the spectral function h(x\,X2) on 
the parameters x\ and X2 is shown in Figs. 15161 As we see, 
the joint spectral function displays the various maxima in 
X\, X2 place that indicate the range of the results of the 
effective three-photon splitting and perfectly confirm the 
results exhibited in Fig. [5] Thus, the above results show 
that the efficient generation of three-photon subharmonic 
mode is possible in this scheme. The parametric gain of 
the three-photon down-conversion in the spectral range 
of the maxima of the function h(xi,%2) is the result of 
the most simultaneous phase-matching of two cascading 
second-order processes. 

In the end of the section we shortly discuss QPM 
for generation of four-photon states in cascaded para- 
metric oscillator. This device consists of three intra- 
cavity modes. The pump mode driven by an external 
coherent driving field at the frequency wq is converted 
into the pair of modes W2, where wq = W2 + VJ2 = 
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FIG. 5: The detailed plot of the function h(xi, X2)in the range 
of double phase-matching peaks: 1 < xi < 2 and 1 < X2 < 2, 
and the parameters N = 22, M = 8. There are four peaks: 
two positive , which are correspond to x\ = 1.5, X2 = 1.5, 
x\ — 1.64, X2 = 1.64, and two negative at x± = 1.5, X2 = 1.64 
and x\ — 1.64, X2 = 1.5. 
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FIG. 6: Plot of the function h(x\,X2) in the ranges 4.2 < 
x\ < 5 and 1 < X2 < 2, and the parameters N — 22, M = 9. 
The heights of the peaks are three times lower than for the 
peaks displayed on Fig. [5] in an agreement with the results 
presented on Fig. [2] 

ujq/2 + u>o/2. Then, each of the mode W2 is transformed 
into second pair of modes w\ , where W2 = roi + w\ = 
n7 o/4 + n7o/4. These processes are characterized by the 
wave-vector mismatches: Afci = k(zuo) — 2k(vD2) and 
Afc 2 = k{vj2) — 2k{w\), respectively. To realize QPM, 
we use the typical value of four-photon down-conversion 
pump wave length X p = 0, 39/im. In LiTa03 phase- 
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FIG. 7: Plot of the function Ym,n{x) for the parameters for si- 
multaneous QPM of four-photon down conversion. The num- 
ber of domains is iV = 32 and that of the blocks is M = 13. 
The insert is detailed the structure of twin-peaks. 



reversed crystal mismatch vectors are Aki = 1,56/mi -1 
Afci = — 1,312/xto -1 , these parameters are taken from 
the experimental results for fourth-harmonic generation 
[l8j . We found the parameters of the phase-reversed 
structure for simultaneous realization of the four-photon 
down-conversion. QPM of mentioned process is realized 
and perfectly matched with the experimental results for 
parameters I — 2, 2/im, M — 13 and N — 32 (see, Fig. [7]). 



IV. CONCLUSION 

In conclusion, we have analyzed phase-reversed QPM 
structures in supcrlattice as the ensemble of domains and 
blocks of second-order nonlinear layers. In the result, the 
effective nonlinear coefficient G(Ak) of three-wave inter- 
action has been obtained as the product of the standard 
phase-matching spectral function for a single segment 
and the assemble function. Generally used phenomeno- 
logical method for description of multilayered structures 
required huge number of domains and blocks to provide 
the periodicity of the function. Our investigations shows 
that the detailed method we have used for description of 
the systems qualitatively matches with phenomenological 
results for the big number of domains N and blocks M, 
moreover, it is applicable for the finite number of layers. 
The results shows that we can QPM various multiphoton 
processes managing number of domains and blocks. We 
have demonstrated that double phase-matched interac- 
tion on the base of phase-reversed QPM structures in su- 
perlattice becomes possible in a wide region of the optical 
wavelengths, particularly, for production of muliphoton 
states. 
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